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ABSTRACT 

We present spectra and high resolution images of the hard X-ray feature along the southern 
edge of the supernova remnant IC443. Data from the Chandra X-ray Observatory reveal 
a comet-shaped nebula of hard emission, which contains a softer point source at its apex. 
We also present A20 cm, A6 cm, and A3. 5 cm images from the Very Large Array that clearly 
show the cometary nebula. Based on the radio and X-ray morphology and spectrum, and 
the radio polarization properties, we argue that this object is a synchrotron nebula powered 
by the compact source that is physically associated with IC443. The spectrum of the soft 
point source is adequately but not uniquely fit by a black body model (kT=0.71±0.08keV, 
L=(6.5 ± 0.9) x 10 31 ergs _1 ). The cometary morphology of the nebula is the result of the 
supersonic motion of the neutron star (Vns — 250±50 kms -1 ), which causes the relativistic 
wind of the pulsar to terminate in a bow shock and trail behind as a synchrotron tail. This 
velocity is consistent with an age of 30,000 years for the SNR and its associated neutron star. 
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Subject headings: ISM: supernova remnants — ISM: individual object: IC 443 — stars: neutron 
-stars: pulsars: individual object: CXOU J061705.3+222127 



Introduction 



The mixed-morphology Galactic supernova remnant IC443 (l,b = 189.1°, +3.0°) has been the subject 



of extensive studies at all wavelength bands (Fesen 1984; Braun fc Strom 1986; van Dishoeck, Janscn fc 



Phillips 1993; Asaoka fc Aschenbach 1994). IC443 is especially well known for its clear interaction with 



surrounding molecular clouds (Burton et al. 1990; Bocchino & Bykov(2000)). With its large variety of 



shocked molecular species detected along its SE edge, IC 443 has become a standard laboratory for studying 
shock chemistry (van Dishoeck et al. 1993 and references therein). 

IC443 is also of interest because it is coincident with the unidentified EGRET source 2EG J0618+2234 
(Sturner & Dermer 1995; Esposito et al. 1996| ), and thus has stimulated a good deal of theoretical work 
aiming to explain the production of GeV 7-rays by shell SNRs (Sturner, Dernier & Mattox 1996; Sturner et 



al. 1997| ; |Gaisscr, Protheroe fc Stanev 1998| paring et al. 1999| ; |Bykov et al. 2000[ ) 



In the X-ray band, Petre et al. (1988) found that the bulk of the emission was thermal in origin with 
a temperature ~10 K, typical of middle-aged SNRs. However, Wang et al. (1992) found evidence in 
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the Ginga satellite data for a hard X-ray component in the spectrum of IC 443 extending up to 20 keV. 
Follow-up observations with the ASCA satellite revealed that the bulk of this hard X-ray emission came 
from a single, unresolved feature located along the edge of the radio shell where the SNR/molecular cloud 



interaction was strongest ( Keohanc et al. 1997 , hereafter K97). Furthermore, the hard X-ray feature is 



positionally coincident with a region where the radio spectral index is considerably flatter than the SNR as 



a whole (i.e a < 0.24 vs a = 0.42, where F v oc v~ a ) (Green 1986; Kovalenko, Pynzar fc Udal'Tsov 1994) 



K97 considered different explanations for the origin of this hard X-ray emission, and favored a model 
in which synchrotron emission was being produced by TeV electrons in a region of enhanced particle 
acceleration, resulting from the SNR/molecular cloud interaction. Recent BeppoSAX data analyzed by 
Bocchino & Bykov (2000) also support the K97 interaction model. An equally viable model, favored by 
Chevalier (1999, hereafter C99), posits that the hard X-ray feature is synchrotron emission powered by 
an energetic neutron star. To distinguish between these two competing hypotheses we have undertaken 
high resolution X-ray and radio observations toward the hard X-ray feature. These new Chandra and VLA 
observations strongly favor C99's neutron star hypothesis. 



2. Observations and Analysis 

2.1. Chandra X-Ray Observatory 

The Chandra X-Ray Observatory performed a short (10 ks) observation of IC443 on April 10, 2000 
during the first cycle of guest observations (A01). The hard X-ray feature was centered on the 13 chip of the 
Advanced CCD Imaging Spectrometer (ACIS). To produce Chandra event files, we followed the standard 
CIAO procedures outlined in Version 1.3 of the CIAO Beginner's Guide for CIAO Release 1.1, using the 
acisD2000-01gainN0001 . f its gain file. From the events file, we extracted high resolution images in both 
hard (E > 2.1 keV) and soft (E < 2.1 keV) spectral bands. Smoothed versions of these images are presented 
in Figure]!]. 

The hard X-ray source is resolved by this high resolution (~1") Chandra observation (Figure [l]) 
as a nebular region of diffuse emission with a cometary tail. Within the nebula there is an unresolved 
point source that we have designated as CXOU 3061105. 3+222127. The nebula of hard emission exhibits 
bow-shock morphology with a width of ~ 35" and a minimum standoff distance of r s ~ 8.5". The X-ray 
point source lies at the apex of the nebula, located at (epoch 2000) a = 06' , 17 m 05.31 s , S = 22°21'27".3, with 
errors of ±2" in each coordinate. It is especially visible in the soft (E < 2.1 keV) X-ray band (Figure]!]). 

Likewise, in order to compare the X-ray luminosity of the cometary nebula with statistical studies of 
pulsar wind nebulae (PWNe) (e.g. Seward & Wang 1988, see §J|), we integrated the power-law model of 
K97 over the Einstein band (0.2-4 keV) - thus deriving L x ~ 5 x 10 33 d 2 (1.5 kpc) -2 erg s^ 1 . Inherent in 
this h x calculation is the assumption that the cometary nebula is dominated in soft X-rays by the same 
power-law exibited at hard energies with ASCA. The Chandra flux (over the same range) yields consistent 
results, albeit with larger uncertainties. 
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Fig. 1. — X-ray intensity images of the hard X-ray feature on the southern edge of IC 443. The image on the 
left is in the soft (E < 2.1 keV) photon energy band, scaled between to 24.6 counts per beam. The image 
on the right is in the hard (E > 2.1 keV) photon energy band, scaled between and 8.5 counts per beam. 
Both of these Chandra images have been smoothed with a 2" beam. The cometary nebula is most apparent 
in the hard band image (right), while the soft point source stands out well against the surrounding nebula 
in the soft band image (left). 



2.2. Very Large Array 

All radio observations were made on August 26 and December 31, 1997 with the Very Large Array 
(VLA0) in the C and D arrays, respectively. A log of the observations is summarized in Table |l|. The data 
acquisition and calibration were standard, using J0632+103 as a phase calibrator and 3C 138 as both the 
flux density and polarization angle calibrator. 

The cometary and bow-shock morphologies of the extended source are even more pronounced at radio 
wavelengths, extending some 2' to the northeast of the soft X-ray source. Although we display only the 
A3. 5 cm image in Figure similar structure is visible at A6 cm and A20 cm. The peak in the radio emission 
lies about 6" to the northeast of the compact X-ray source. There is no radio point source (> 2 mJy) 
coincident with the X-ray point source. The distribution of linearly polarized emission is also shown in 
Figure ||. The direction of the E- vectors are generally parallel to the shock normal, as would be expected for 
evolved SNRs with circumferential magnetic fields. The character of the magnetic field evidently changes 
substantially in the area dominated by the X-ray emission. Here the magnetic field direction appears 
to wrap around the "head" of the hard nebula, tracing a bow-shock morphology. The average degree of 
polarization in this region is 8%. This is likely a lower limit since, owing to the compact size of the nebula, 
there is likely some beam depolarization occurring. To the northeast, where the field is more ordered, 
the percentage polarization is uncharacteristically strong for an evolved SNR, exceeding 25% in several 
locations. 
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Fig. 2. — Total intensity image of the cometary nebula in IC 443 (left) at a frequency of 8.46 GHz. Contours 
are in steps of 0.5 mJy beam -1 . The size of the synthesized beam is shown in the lower left corner. The 
same total intensity image with vectors indicating the magnitude and direction of the polarized intensity 
(right). The scale is 50 fiJy beam -1 arcsec -1 . The cross marks the position of CXOU J061705. 3+222127. 

The integrated flux density over the entire length of the cometary nebula is given for three frequencies 
in Table |l} An approximate value of 230±200 mJy was also obtained from the 327 MHz image presented 
by Claussen et al. (1997). The increased error in estimating the flux density at low frequencies is due to 
the uncertainty in subtracting out the background shell emission from the SNR. Within the errors, the 
radio spectrum is well represented by a flat spectral index an ~ 0.0 and a mean flux density of 206 mJy. 
In contrast, K97 measured an X-ray spectral index of ax ~1.3±0.2, which implies a break in the spectrum 
near 4xl0 13 Hz. The radio luminosity of the cometary nebula, determined by integrating the spectrum 
from 10 MHz to 100 GHz, is L fl =5.5xl0 31 (d/1.5kpc) 2 ergs -1 , or less than 1% of L R for IC443 as a whole. 



Table 1. Summary of VLA Observations of the Cometary Nebula 



Frequency 


Time 


Beam 




(GHz) 


(min) 


(" x ") 


(mJy) 


8.46 


13 


8.6 x 7.6 


195 ± 8 


4.86 


36 


5.0 x 4.8 


173 ± 11 


1.46 


53 


15.5 x 14.5 


229 ± 34 



Note. — The columns are (left to right), (1) Observing frequency, (2) total time on source, (3) angular resolution, and (4) 
integrated flux density of the cometary nebula. 



- 5 - 



CXOCO J061705. 3 + 222127 with Black Body Model 
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Fig. 3. — The background subtracted spectrum of the soft point source, fit with a black body emission model 
(kT=0.7±0.1 keV and A=0.02 ± 0.01(d/1.5 kpc) 2 km 2 ), assuming a previously measured column density. 
Channels 5, 10 and 15 correspond to photon energies of 1.5, 2.9, and 4.2 keV respectively. 



2.3. Spectra 

In addition to the spectral properties of the hard nebula, a spectrum of CXOUJ061705. 3+222127 was 
extracted from this short Chandra observation. The surrounding hard nebula was subtracted as background 
(i.e. an annulus with an inner and outer radius of 4.3" and 23.6" respectively). These spectra were binned 
by a factor of 256 (PHA channels), and the softest three ACIST channels were ignored. As discussed above, 
the point source is significantly softer than the surrounding hard nebula, implying a likely thermal origin 
for the emission. For this reason, we modeled the spectrum as black body radiation, though this fit was not 
statistically unique due to large uncertainties in the small data set. 

The best-fit (xl ~ 1-1) black body temperature and luminosity are kT=0.71±0.08 keV and 
L=(6.5 ± 0.9) x lO^ergs" 1 (i.e. A=0.025 ± 0.003(dl.5 kpc) 2 km 2 ). The 1-5 keV X-ray flux of the point 
source is ~ 2 x 10 -13 ergs -1 cm~ 2 . This fits contained two free parameters by assuming K97's best-fit 
column density (NH = 1.3xl0 21 cm -2 ). 



Discussion 



The soft X-ray point source in IC 443 is one of a growing number of young neutron star candidates 



associated with supernova remnants ( Hclfand 1998 ). These neutron stars have a diverse range of properties, 



suggesting that they do not originate from a homogeneous population (Chakrabarty et al. 200C). The 
long-period (8-10 s) and claimed high-field (10 14 -10 15 G) of anomalous X-ray pulsars (AXP) and soft 
gamma-ray repeaters (SGR), differ from the canonical young radio pulsars with periods P~15-500 ms and 
field strengths of B~10 12 -10 13 G. The radio quiet neutron stars (RQNs), with their thermal X-ray emission 



and lack of radio pulsations, may be altogether different from either of these two classes of objects ( Pavlov 



et al. 20001) 



In the absence of detectable pulsations, canonical radio pulsars can be distinguished from AXPs, 
SGRs, and RQNs by the presence of an extended synchrotron nebula, powered by the energetic relativistic 
wind ( |Gacnslcr, Bock fc Stappers 2000 ). It is argued that the cometary nebula around the X-ray point 
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source in IC 443 is most likely such a PWN. This nebula has all the expected observational characteristics 



(C99; |Gaensler et al. 200C): flat spectrum radio emission (or = 0.1 to 0.3), a steep X-ray spectrum {ax 



1.0 to 1.5) and a high degree of linear polarization (> 5%). Morphologically it closely resembles PWN 



detected toward the SNRs W44 and G 5.4-1.2 ( [Frail fe Kulkarni 199J ; [Frail et al. 199(i| ). both of which 



contain active pulsars. It has been argued (Frail & Kulkarni 1991; Frail et al. 1996) that these nebulae have 



developed bow shocks as a result of ram pressure confinement of the pulsar wind, due to the high space 
velocity of their pulsar through the surrounding medium. We assert that the most tenable explanation for 
the X-ray point source in IC 443 is that it is a young pulsar which has traveled (ballistically) from its birth 
place to its present location, near the edge of the decelerating SNR. 

Accepting this hypothesis we proceed with inferring the salient physical properties of the neutron star 
and corresponding PWN. Constraints can be placed on the velocity of the neutron star by a variety of 
methods. The offset of the object from its origin gives a measure of its transverse velocity Vns- Measuring 
this offset in the case of IC 443 is made difficult owing to the uncertain location of the explosion center 
( Asaoka k, Aschenbach 1994 ; Claussen et al. 1997] ) . We note that the neutron star has traveled a fraction 



/3 ~0.9 of this projected distance (approximated from the distance between the soft point source and 
the shock), giving Vns = /3r s /i s . Likewise, the shock velocity of the SNR can be similarly expressed as 
V s = c T s /t s , where c is a constant equal to 2/5 for a remnant in the Sedov stage, while C99 argues for a 
value of c = 3/10 for IC443 based off of the parameters of his model. Combining both equations for Vns 
and V s together, and using the value of V s = 100 kms" 1 favored by C99, gives Vns = V s /3c o _1 =225-300 
km s -1 , a result that is independent of distance and only weakly dependent on the evolutionary state of the 
SNR (i.e. Co). 

It should be noted that the tail of the hard nebula does not point towards the geometric center of the 
SNR. However, the blast center and geometric center of an SNR can be quite different in the presence of 
large-scale density gradients ( Dohm-Palmcr fc Jones 1996] : Hnatyk fc Pctruk 199E ). Such a density gradient, 



combined with the potential of a cross-wind, could possibly account for the apparent discrepancy. 

We are able to make a second approximation of Vns with respect to the local medium by means of 
the nebula's bow shock morphology. K97 fit the SNR with an X-ray temperature kT ~ 1 keV, implying a 
sound speed Ci oca i ~100 kms" 1 . An approximation of the bow shock angle of 90°-120° gives a minimum 
Mach number M~l.l-1.5 through a standard geometric calculation. A fundamental determination of the 
speed of sound in this particular region of the SNR gives velocity with respect to the local medium of 
VNSiocai=T50Mi.5 kT^v kms -1 . Accounting for radial expansion of local plasma (<100 kms" 1 , this 
velocity could be as high as 250 kms -1 . Both methods of determining Vns are m good agreement with a 
model by C99, indicating a remnant age on the order of 30,000 years (i.e. T30 ~1). 

The energy source for the synchrotron nebula must ultimately be derived from the particles and field 
generated by the compact central neutron star and therefore some fraction of the its spindown power E is 
powering this emission. Empirical relations have been derived for known PWN between E and the observed 
X-ray and radio luminosities (L x and L#) that enable us to estimate the rotation period P and the magnetic 
field B of the compact object. 

For the value of h x derived in §|| we obtain E~2xl0 36 ergs -1 (Seward and Wang 1988). Frail & 
Scharringhausen (1997) and Gaensler et al. (2000) note that E~10 4 L# for known radio PWN, which yields 
a value of E~6xl0 35 , so in the discussion that follows we adopt E=10 36 ergss" 1 (i.e. Fj36=1). Equating 
the pressure of the wind E/47rcr^, to the ram pressure due to the neutron star's motion through the 
surrounding gas pVp SR , we derive a density of ~0.1 cm" 3 for our adopted values of E and Vns- This 
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density is more typical of that expected for the hot X-ray interior of the SNR, and is orders of magnitude 
below that expected in the dense molecular ring against which the NS and its nebula are projected against 
(see van Dishoeck et al. 1993). 

The spindown luminosity of a pulsar is related to its current period P and dipolar surface held B by 
EocB 2 P~ 4 , while the characteristic age of the pulsar is expressed as t c ocB~ 2 P 2 . It is straightforward to 
show that for a pulsar born spinning rapidly (i.e. P <<P), which loses energy predominately via magnetic 
dipole radiation, P=145ms (t 30 E36)~ 1/2 and B=3.3xl0 12 G t^E"^ 2 . Thus, with our nominal choice of 
parameters, the period and magnetic field closely match those of other young pulsars (e.g. PSR B1757— 24; 



Frail fc Kulkarni 1991) 



In the light of this result, it is worth re-examining the nature of the unidentified EGRET source 
2EG J0618+2234 (l,b=189.13°,+3.19°) flSturncr fc Dormer 1995] |Esposito et al. 1996|) toward IC443. Nel 



et al. (1996) note that all known gamma-ray pulsars have a large (greater than 0.5) ratio of E^Jdkpc ■ 
For our derived E ~ 10 36 ergs _1 and d^l.5kpc, this ratio is over 400, similar to that of the gamma-ray 
emitting PSR B1706— 44, strongly suggesting that the unidentified EGRET emission originates from the 
neutron star or the PWN. The difficulty with this interpretation is that our proposed neutron star and 
PWN lie several arcminutes outside the formal 95% error radius of the EGRET source. However, since the 



position of EGRET sources in the Galactic plane are subject to systematic error (e.g. Hunter et al. 1997 ), 
it is entirely possible that this object may be the source of this gamma-ray emission. 



4. Conclusion 

In summary, the hard X-ray feature in IC 443 is best interpreted as a wind nebula powered by a 
young neutron star. The PWN interpretation is suggested first on morphological grounds by the Chandra 
and VLA images, which show a comet-shaped nebula with a soft X-ray point source at its apex. The 
measurement of significant polarization and a flat radio spectrum further strengthens the argument that 
this is synchrotron emission from a PWN. The X-ray spectrum of the point source is appropriately, but not 
uniquely fit by a black body model, suggesting that the emission is thermal in origin. The inferred physical 
properties of the nebula and point source (i.e. r, Vns> E, P, B) support our hypothesis that IC443 and the 
pulsar were produced approximately 30,000 years ago in a core collapse supernova event. 

Despite the evidence in favor of a real physical association between IC 443 and the compact source 
within the cometary nebula, there are a few puzzles that remain to be explained. As mentioned previously, 
one might expect the "tail" of the cometary nebula to point back to the origin of the SNR and the pulsar. 
However, the geometric center of IC443 lies almost due north of the compact object, and not to the 
northeast as suggested by the direction of the cometary tail. Given the complex kinematics and distribution 



of the molecular and neutral gas in the vicinity of IC443 (Giovanelli & Haynes 1979; van Dishoeck et al. 
1993) it may be overly naive to ascribe the blast center to the geometric center. A final problem concerns 
the apparent absence of X-ray and radio pulsations from the IC 443 point source. K97 carried out a pulse 
search using ROSAT and ASCA data, while Kaspi et al. (1996) have searched for radio pulsations. In our 
view, these null results present no immediate problem for our young neutron star hypothesis since pulsed 
X-rays from Vela-like pulsars have proven difficult to detect due to their soft spectra and contamination 



from the surrounding nebular emission (e.g. Becker & Trucmpcr 1997). Furthermore, it is entirely possible 



that the radio beam may not intersect our line of sight (Lorimcr, Lync & Camilo 199 
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